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Abstract 



Quantum Chromodynamics (QCD), the gauge field theory of the Strong Interaction, has specific 
features, asymptotic freedom and confinement, which determine the behaviour of quarks and gluons 
in particle reactions at high and at low energy scales. QCD predicts that the strong coupling 
strength as decreases with increasing energy or momentum transfer, and vanishes at asymptotically 
high energies. In this review, the history and the status of experimental tests of asymptotic 
freedom are summarised. The world summary of measurements of Os is updated, leading to an 
unambiguous verification of the running of Qg and of asymptotic freedom, in excellent agreement 
with the predictions of QCD. Averaging a set of measurements balanced between different particle 
processes and the available energy range, results in a new and improved world average of as(M2o) = 
0.1189 ±0.0010 . 



MPP-2006-54 
June 2006 



1 



Contents 



1 Introduction 



2 Historical Development 



2.1 A New Force! 



2.2 



Quarksl 



2.3 Chromo-Dvnamicd 



2.4 Asymptotic Freedom 



2.5 Scaling Violations! 



2.6 Quark and gluon iets in high energy collision^ 



2.7 Status in the late 1980's! 



3 Theoretical Basis! 



3.1 Renormalization 



3.2 and its energy dependence! 



3.3 Asymptotic freedom and confinement 



3.4 The running coupling 



3.5 Relative size of finite order approximation^ 



3.6 Quark masses and thresholds! 



3.7 Perturbatiye predictions of physical quantities! 



3.8 Renormalization scale dependence, 



3.9 Nonperturbatiye QCD methods! 



4 Tests of Asymptotic freedom in e^e annihilations! 



4.1 Jet production rates! 



4.2 Evidence for the gluon self coupling 



4.3 Determination of the QCD group constant^ 



4.4 The running from 4-iet event production 



5 QCD tests in deep inelastic lepton-nucleon scattering! 



5.1 Basic introduction to structure functions! 



5.2 Scaling violations of structure functions! 



5.3 The running a= from iet production in "pis! 



6 Summary of a„ measurements! 



6.1 Previous developments and current status! 



6.2 Update: new results since 2004 



6.3 Qc. summary 



6.4 A new world average of a^iM^.o) 



7 Summary and Outlook 



1 Introduction 



The Strong Interaction, binding quarks and gluons inside hadrons, is the strongest of the four funda- 
mental forces in nature which we know today. It is about 100 times stronger than the electromagnetic 
force, a factor of 10^^ stronger than the Weak Interaction, and a stunning factor of lO^'' stronger^ than 
the Gravitational Force, calculated for two quarks at a distance of order 1 fm. 

In spite of these differences, it is the Gravitational and the Electromagnetic Forces which seem to 
play the dominating role in the universe which we experience in the macroscopic world. The reason for 
this wondrous fact is that the Strong Force, as well as the Weak Interaction, only acts at subatomic 
distances. Taking the existence of atomic nuclei, their composition and their masses for granted, 
restricting ourselves to energies and temperatures below 1 MeV and, respectively, 10^ Kelvin, the 
Strong Interaction is disguised in the description of the world. 

At subatomic distances, however, this picture changes completely. The Strong Force not only 
determines the binding of quarks and gluons inside hadrons, it also determines the cohesion of protons 
and neutrons inside atomic nuclei. Hadrons like protons and neutrons are responsible for more than 
99% of the mass of all visible matter in our universe, and those masses are mainly generated by the 
strong binding of quarks inside hadrons, rather than by the (generally small) masses of the quarks 
themselves. 

The restriction of the Strong Force to subatomic distances is a consequence of two characteristic 
features: they are called "Confinement" and "Asymptotic Freedom". 

Confinement is a necessary requirement to explain the fact that no isolated quarks have ever been 
observed in any experiment, although symmetry arguments and scattering experiments in the 1960's 
established quarks, with -1/3 or +2/3 of electrical charge units and a newly introduced quantum 
property called "colour charge", as the basic constituents of hadrons. Confinement determines that at 
large distances, or - equivalently - at low momentum or energy transfer in elementary particle reactions, 
the Strong Force prevents the existence of free quarks: Trying to separate two quarks from each other, 
for instance in high energy scattering reactions, apparently results in an increase of the force field's 
energy at large distances, such that new quarks are created out of the vacuum - the initial quarks 
"dress" up with other quarks to build hadrons. These hadrons exhibit no net colour charge to the 
outside, such that they appear as elementary entities rather than the quarks themselves^. 

The term "Asymptotic Freedom" is used to describe the behaviour of quarks at high energy or 
momentum transfers, or - equivalently - at small distances. Also this feature is based on experimental 
observations: In high energy scattering processes between leptons (e.g. electrons or neutrinos) with 
protons or neutrons, the dynamics reveal that scattering occurs at pointlike and massless constituents, 
the quarks, rather than at a homogenuous object with the size of a proton. Apparently, at sufficiently 
high momentum transfers, quarks behave like free or weakly bound particles. Also, quarks knocked out 
of a hadron, in a high energy scattering process, were never observed as free particles. Instead, they 
emerge as dressed-up hadrons or bundles (jets) of hadrons escaping from the interaction region. 

The fact that the strong interaction becomes "weak" at high energy scales, and vanishes to zero 
at asymptotically high energies, led to the term "Asymptotic Freedom". Any theory describing and 
predicting the dynamics of quarks inside hadrons and in high energy reactions therefore had to satisfy 
and include these two extremes: Confinement, also called "infrared slavery", and Asymptotic Freedom. 

A theoretical description of the Strong Interaction by a consistent quantum field theory, called 
Quantum Chromo dynamics (QCD), was presented in 1972 by Fritzsch and Gell-Mann pQ, which was 
formaly published later in |2] . At the same time it was found by Gross, Politzer and Wilczek that non- 
Abelian gauge theories, such as QCD, exhibit the propoerty of Asymptotic Freedom This discovery 
was honored with the Nobel Price in 2004, 30 years after the original findings. 

"'^40 orders of magnitude correspond to the difference in size of our entire universe to the size of an atomic nucleus. 
^Van-der-Vaals-hke remnant forces bind protons and neutrons in atomic nuclei. 
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Theoretical breakthroughs and discoveries, before being honored by the Nobel Price, must be es- 
tablished by experimental measurements. The price for Asymptotic Freedom therefore also had to 
await experimental verification. In this article, the long but exciting road to experimental evidence for 
Asymptotic Freedom shall be summarised. 

2 Historical Development 

2.1 A New Force 

In the 1930's, protons and Neutrons were recognised as the building blocks of atomic nuclei. A new and 
elementary force, the so-called Strong Force, was introduced phenomenologically, in order to explain 
the binding of these particles inside nuclei. These binding energies were known to be much larger than 
what could be expected from the well-known electromagnetic force. Also, the electromagnetic force was 
known to have infinite range while the Strong Force apparently acted only within nuclear distances of 
the order of femtometer. 

Shortly after these perceptions, Hideki Yukawa argued that the range of the electromagnetic force 
was infinite because the associated exchange particle, the photon, was massless, and proposed that the 
short range of the strong force was due to the exchange of a massive particle which was called "meson" 

From quantum mechanical arguments and from Heisenberg's uncertainty principle Yukawa predicted 
a meson mass^ of about 100 MeV/c^. At that time, no particle with a mass in that range was known. 

Very soon after, in 1937, a new particle was dicovered whose mass was very close to that of Yukawa's 
prediction of the meson - however, this particle in fact was the muon which turned out not to be subject 
to the strong interaction. Finally, in 1947 strongly interacting particles with masses close to Yukawa's 
prediction were found [Sj in emulsion experiments at high altitudes. A particle which met Yukawa's 
requirements was found, and it was called "pion". Although today we know that pions, like protons 
and neutrons, are composite particles made from quarks, and the Strong Interaction is not mediated 
through the exchange of pions, Yukawa's theory stimulated major advances in the understanding and 
description of the strong interaction. 

2.2 Quarks 

In the decades after World War II, the number of strongly interacting particles exploded - which was 
basically due to the development of particle accelerators with steadily increasing energies. Strongly 
interacting particles were then called "hadrons" , which differentiate into mesons (with integer or zero 
spins) and baryons (with spins of 1/2 or 3/2). Gell-Mann and Zweig realised in 1964 that the whole 
spectroscopy of hadrons could be explained by a small number of "quarks" if baryons are made out of 3 
quarks, and mesons out of one quark plus one antiquark P . Those quarks, initially 3 different species, 
then must have 1/3 or 2/3 of the elementary charge unit, and they must have spin 1/2. 

By the end of the 1960's, the static picture of quarks as the constituents of hadrons was confirmed 
through the dynamics observed at high energy electron-proton scattering experiments at the Stanford 
Linear Accelerator (SLAG) 7j, summarised in figure ^ the cross sections - instead of decreasing with 
increasing momentum transfer as expected for elastic scattering of electrons at protons as a whole - 
showed a "scaling" behaviour, as it should occur if the electrons scatter on quasi-free, pointlike and 
nearly massless constituents inside the protons. 

While the quark model was very successful in describing the properties, multitude and dynamic 
behaviour of hadrons, it had severe shortcomings such as an obvious violation of the Pauli-principle 

•^From now on, in this review, a system of units is utilized where the speed of hght and Planck's constant are set to 
unity, c = h = 1, such that energies, momenta and masses are all given in units of eV. 
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Figure 1: The structure function of the proton, describing the electron- proton scattering in units of 
the Mott cross section, i.e. the generahsed Rutherford cross section at high energies, as measured at 
SLAC uj = 4 corresponds to x = 0.25 in figure IT^ 



for hadrons containing two or three identical quarks in relative S-waves and in identical spin states, 
the prediction of the neutral pion's lifetime which was wrong by a factor of nine, and the fact that 
no particles (quarks) with 1/3 or 2/3 of the elementary charge unit could ever be observed at particle 
colliders. 



2.3 Chromo-Dynamics 

To overcome these shortcomings, the idea that the strong force couples to a new quantum poperty, 
called colour charge jHl E] , in analogy to the photon in QED which couples to the electric charges. 

The introduction of the new quantum number, describing 3 different colour quantum states of each 
of the quark species, solved the puzzle of spin-statistics, saved the Pauli-principle and explained the 
missing factor of nine (= 3^) of the pion lifetime. The notion that hadrons consist either of three quarks 
(baryons like the proton) or a quark and an antiquark (mesons), arranged such that the net color charge 
of the hadron would vanish - i.e. "white" as a superposition of the three elementary color states or of 
a particular color and its anti-color - could account for the fact that the strong force is short-ranged. 

Finally, in the early 1970's, a field theory of the strong force. Quantum Chromodynamics (QCD), was 
introduced and presented in canonical form P3I2]- Here, coloured spin-1 particles called "gluons", which 
- in contrast to the case of photons in QED - carry colour charge themselves, couple to the colour charges 
of quarks, and also to coloured gluons themselves. Chromo-Statics turned into Chromo-Dynamics. 



2.4 Asymptotic Freedom 

A few years before, it was demonstrated in model theories that charges may change their effective size 
when they are probed in scattering experiments, at large and at small distances. These changes were 
described through Symanzik's /3-function [TU], see equation 2 below: the effective size of the coupling 
strength is a function of the energy or momentum transfer, Q^. The SLAC data on approximate scaling 
of the proton structure function, c.f. figure ^ and the notion of free quarks inside the proton required 
a negative /3-function. At that time, however, all field theories probed so far had a positive /^-function. 
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Gross, Politzer and Wilczek finally demonstrated in 1973 that Chromo-Dynamics, with coloured 
quarks and gluons, obeying S'U{3)coiour symmetry, generated a negative /3-function i.e. that quarks 
and gluons are asymptotically free. This explained the SLAG data, and at the same time leads - for 
small momentum transfers or at large distances, to an increase of the coupling strength, thus motivating 
confinement. 

Another important consequence from asymptotic freedom is the fact that the strong coupling ag is 
small enough, at sufficiently large momentum transfers, to allow application of perturbation theory in 
order to provide quantitative predictions of physical processes. 

With this discovery Quantum Ghromodynamics started its triumphal procession as being the field 
theory of the strong interaction. It was, however, a long and difficult path before QGD was commonly 
accepted to be exactly that. Many refined calculations, theoretical predictions and experimental verifi- 
cations were ventured. The most important features of QGD, asymptotic freedom and/or, equivalently, 
the existence of colour charged gluons, had to be tested, quantified and "proven" by experiment. The 
size of the strong coupling parameter, a^iQ^), had to be determined and its energy dependence verified 
to be compatible with asymptotic freedom. 

2.5 Scaling Violations 

The violation of approximate scaling of the proton structure function was one of the first experimental 
signatures proposed to test QGD |ir . Deep inelastic electron-proton (or more generally, lepton-nucleon) 
scattering processes basically depend on two kinematical parameters. The usual choice of these param- 
eters are the momentum transfer between the lepton and the struck quark, and the fraction x of the 
proton's momentum which is carried by the struck quark. 

The structure function F2 of the proton basically parametrises the population of quarks with momen- 
tum fraction x, as a function of Q^: F2{x, Q"^). QGD predicts that the region of large x is depopulated 
for higher momentum transfers, while the population at low x should increase with increasing . 

A simple argument for this behaviour can be given in terms of the equivalence of momentum transfer 
and physical resolution: with higher Q^, smaller distances are probed. At smaller distances, an increased 
radiation of soft (low energy) gluons should be resolved. These gluons increase the population of partons 
at small x, and - at the same time - diminish the momentum fraction and the relative contribution of 
quarks at large x. 

On top of it's kinematical origin, scaling violations in x and are modified by the specific energy 
dependence of the strong coupling. Note that a qualitative observation of scaling violations of structure 
functions alone does not yet prove Asymptotic Freedom - such a conclusion requires a large lever arm 
in and a precise study of the functional form of observed scaling violations, as e.g. the logarithmic 
slopes, dF2{x,Q'^)/dlnQ^. 

While scaling violations were already observed in the 1970's, initially it was not possible to conclude 
that asymptotic freedom was actually seen |T21IISI- Instead, unknown higher order QGD terms, the 
effects of target masses and nonperturbative effects ("higher twists") could cause similar effects as those 
predicted by asymptotic freedom. Only with the availability of much higher and the ability to probe 
much smaller values of x, e.g. at the electron-proton collider HERA, stricter conclusions were possible 
in the past 10 to 15 years. These will be presented in more detail in section 5. 

2.6 Quark and gluon jets in high energy collisions 

Another class of predictions, basically governed by the concept of confinement, was the expectation of 
collimated hadron jets in high energy reactions like e^e~ qq Indeed, by 1975, the emergence 

of 2-jet structure was observed when increasing the e^e~ center of mass energy Ecm from 3 to 7.4 GeV 
|16j . confirming the basic ideas of the quark-parton model. 
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The radiation of energetic gluons off high energetic quarks was predicted from QCD, and hence the 
emergence of a third hadron jet was expected in e"'"e~ annihilation at higher energies jT7]. In 1979, 3-jet 
structures were observed at the PETRA e"'"e~ coUider at Ecm ~ 30 GeV [TH], see figure |21 The 3^*^ jet 
could be attributed to the emission of a third parton with zero electric charge and spin 1 [T^ - the 
gluon was discovered and explicitly seen! 



Figure 2: The first 3-jet event, observed by the TASSO experiment at the PETRA e"''e storage ring 

CHI- 



After an almost parameter-free description of hadronic event shapes [201 1^ had served to establish 
gluon radiation beyond any doubts, first determinations of the coupling strength as were performed. 
More detailed studies of the jet structure and the shape of hadronic events provided insights into the 
color structure of the gluon, like for instance the "string effect" |221 which was expected in QCD due 
to the higher color charge of the gluon [23 • 



Towards the end of the 1980's, after successful exploitation of the experimental programs of the e+e~ 
storage rings PETRA at DESY and PEP at SLAG, the proton-antiproton collider at CERN and many 
deep inelastic lepton-nucleon experiments at CERN, at Fermilab and at SLAC, "t/ie experimental sup- 
port for QCD is quite solid and quantitative" [21], however convincing proofs of the key features of 
QCD, of the gluon-selfcoupling and/or of asymptotic freedom, were still not available. As an example, 
the summary of measurements of ag at different energy scales in 1989, as reproduced in figure El 



A tracks • ■ 
A.I GeV 



■ •■..5 tracks 



A.3 GeV 




2.7 Status in the late 1980's 
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was compatible with the QCD expectation of the running of as, but did not yet allow to draw more 
concrete conclusions. 
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Figure 3: Summary of measurements of as in 1989 [24. Shown are results from various experiments in 
deep inelastic lepton-nucleon scattering as well as combined results from e^e~ collisions, together with 
the QCD expectation of a running as for different values of (see section 3). 

Foundations were layed, however, to prepare for more quantitative tests with the upcoming higher 
energy colliders, like the e^e~ colliders LEP at CERN and SLC at SLAC, the Tevatron pp collider at 
Fermilab and the HERA electron-proton collider at DESY, who all started operation in the time from 



• From the summary of ag measurements shown in figure El using the QCD prediction of a running 
coupling, the value of asi^Q"^ = M^o), at the energy scale of the Z° boson mass, M^o 90 GeV, 
was predicted to be as{Mzo) = 0.11 ± 0.01. Experiments at LEP and SLC were determined to 
scrutinise this prediction with high accuracy 

• Deep inelastic scattering results, extending to much higher values of and much lower values of 
X, should proove scaling violations of nucleon structure functions. 

• Significant progress in jet physics allowed to prepare for direct tests of asymptotic freedom, through 
the energy dependence of jet production rates [23 122] , and of the gluon-selfcoupling, through spin- 
correlations in 4-jet hadronic final states in e"'"e~ annihilation |27j . 

• Precise determinations of as at low and at high energy scales, in e"'"e~ annihilation, at hadron 
colliders and in deep inelastic lepton-nucleon collisions were expected to emerge, with the prospect 
of proving the energy dependence of as, and thus, asymptotic freedom. 



1987 to 1991: 
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• Theoretical developments like precise calculations in next-to- leading (NLO) and next-next-to- 
leading order (NNLO) pertubation theory, advanced methods of handling theoretical uncertainties, 
and improvement in the understanding of the nonperturbative hadronisation process, in terms of 
advanced Monte Carlo models, had started and were expected to be confronted with future high 
precision experimental data. 

After a short summary of the theoretical basics of QCD in section 3, the developments of experimen- 
tal test of QCD, and in particular of asymptotic freedom, following the status as summarised above, 
will be presented in more detail in sections 4 to 6. 

3 Theoretical Basis 

The concepts of QCD are portrayed in many text books and articles, se e.g. |2H1 1211 IHOl IHH and 
in particular, the article of G.M.Prosperi, M. Raciti and C.Simolo in this issue of PPNP [S3]- In the 
following, a brief summary of the basics of perturbative QCD and of the energy dependent strong 
coupling parameter, a^, is given. 

3.1 Renormalization 

In quantum field theories like Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED), 
physical quantities TZ can be expressed by a perturbation series in powers of the coupling parameter 
tts or a, respectively. If these couplings are sufficiently small, i.e. if ag <C 1, the series may converge 
sufficiently quickly such that it provides a realistic prediction of TZ even if only a limited number of 
perturbative orders will be known. 

In QCD, examples of such quantities are cross sections, decay rates, jet production rates or hadronic 
event shapes. Consider TZ being dimensionless and depending on as and on a single energy scale Q. 
This scale shall be larger than any other relevant, dimensional parameter such as quark masses. In the 
following, these masses are therefore set to zero. 

When calculating 7^ as a perturbation series of a pointlike field theory in as, ultraviolet divergencies 
occur. These divergencies are removed by the "renormalisation" of a small set of physical parameters. 
Fixing these parameters at a given scale and absorbing this way the ultraviolet divergencies, introduces 
a second but artificial momentum or energy scale fi. As a consequence of this procedure, TZ and 
become functions of the renormalization scale /x. Since TZ is dimensionless, we assume that it only 
depends on the ratio Q^/fJ,"^ and on the renormalized coupling as(yU^): 

7^ = n{Q'^/iJ, as); «s = as(/i^)- 

Because the choice of fi is arbitrary, however, the actual value of the experimental observable TZ 
cannot depend on fi, so that 

,^_i,K,<^v,^„,).(,^|,,,^|J^)^ao, (1, 

where the derivative is multiplied with in order to keep the expression dimensionless. Equation ^ 
implies that any explicite dependence of 7?. on /i must be cancelled by an appropriate /i-dependence of 
tts to all orders. It would therefore be natural to identify the renormalization scale with the physical 
energy scale of the process, /x^ = Q^, eliminating the uncomfortable presence of a second and unspecified 
scale. In this case, as transforms to the "running coupling constant" as(Q^), and the energy dependence 
of TZ enters only through the energy dependence of as{Q^). 

Any residual /i-dependence is a measure of the quality of a given calculation in finite perturbative 
order. 



9 



3.2 as and its energy dependence 

While QCD does not predict the absolute size of as, its energy dependence is precisely determined. 
If as(/i^) is measured at a given scale, QCD definitely predicts its size at any other energy scale Q"^ 
through the renormalization group equation 

The perturbative expansion of the (3 function is calculated to complete 4- loop approximation |34j : 

/3(a,(Q2)) = -p,a',{Q^) - Aaf (Q') - P2atiQ') - PsaUQ') + Oiaf) , (3) 



where 

/92 



33 - 2Nf 
153 - 19Nf 

77139 - 15099A^/ + 325A^| 
34567r3 



29243 - 6946.3Arf + 405.089iV? + 1.49931Ar3 
* 256.^ ' '-• W 

and Nf is the number of active quark flavours at the energy scale Q. The numerical constants in 
equation m are functions of the group constants Ca = N and Cp = {N'^ — 1)/2A^, for theories exhibiting 
SU{N) symmetry; for QCD and SU{3), Ca = 3 and Cp = 4/3. Pq and /?i are independent of the 
renormalization scheme, while all higher order (3 coefficients are scheme dependent. 

3.3 Asymptotic freedom and confinement 

A solution of equation |2l in 1-loop approximation, i.e. neglecting f3i and higher order terms, is 

l + as(/x2)/3oln^ 

Apart from giving a relation between the values of ag at two different energy scales and /x^, equationEl 
also demonstrates the property of asymptotic freedom: if becomes large and (3q is positive, i.e. if 
Nf < 17, as(Q^) will asymptotically decrease to zero. 

Likewise, equation El indicates that q;s(Q^) grows to large values and, in this perturbative form, 
actually diverges to infinity at small Q^: for instance, with as(/i^ = M|o) = 0.12 and for typical values 
of Nf = 2 ... 5, as{Q^) exceeds unity for < 0(100 MeV ... 1 GeV). Clearly, this is the region where 
perturbative expansions in as are not meaningful anymore, and we may regard energy scales below the 
order of 1 GeV as the nonperturbative region where confinement sets in, and where equations El and El 
cannot be applied. 

Including f3i and higher order terms, similar but more complicated relations for as((5^), as a function 
of as(/i^) and of In % as in equation emerge. They can be solved numerically, such that for a given 
value of as(/x^), choosing a suitable reference scale like the mass of the Z° boson, /i = M^o, as(Q^) can 
be accurately determined at any energy scale > I GeV^. 

If we set 

A^: 
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a dimensional parameter A is introduced such that equation El transforms into 

= Alu(QVA=) ■ 

Hence, the A parameter is technically identical to the energy scale Q where OsiQ"^) diverges to infinity, 
as{Q^) oo for ^ a^. To give a numerical example, A ^ 0.1 GeV for ^^(Mzo = 91.2 GeV) = 0.12 
and Nf = 5. 

The parametrization of the running coupling as(Q^) with A instead of Q;s(yU^) has become a common 
standard, see e.g. [21], and will also be adopted here. While being a convenient and well-used choice, 
however, this parametrization has several subtleties: 

First, requiring that as(Q^) must be continuous when crossing a quark threshold^, A actually depends 
on the number of active quark fiavours. Secondly, A depends on the renormalization scheme, see e.g. 
reference [22] • In this review, the so-called "modified minimal subtraction scheme" (MS) [32] will be 
adopted, which also has become a common standard [21] ■ A will therefore be labelled A^^'' to indicate 
these peculiarities. 
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Figure 4: (a) The running of as{Q), according to equation [71 in 4- loop approximation, for different 
values of Ajj^; (b) same as full line in (a), but as function of l/log(Q/GeV) to demonstrate asymptotic 
freedom, i.e. as(Q^) — > for Q ^ cc. 



3.4 The running coupling 

In complete 4-loop approximation and using the A-parametrization, the running coupling is thus given 
by 



^Strictly speaking, physical observables TZ rather than must be continuous, which may lead to small discontinuities 
in UsiQ^) at quark thresholds in finite order perturbation theory; see section 3.7. 
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+ 



+ 



^-ln'L + -ln^L + 2lnL-- 



(7) 



where L = Q'^/Aj^. The first hne of equation [7| includes the 1- and the 2-loop coefficients, the second 
hne is the 3-loop and the third line is the 4-Ioop correction, respectively. 
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Figure 5: (a) The running of as{Q), according to equation [71 in 1-, 2- and 3-loop approximation, for 
Nf = 5 and the same value of = 0.22 GeV. The 4-loop prediction is indistinguishable from the 
3-loop curve, (b) Fractional difference between the 4-loop and the 1-, 2- and 3-loop presentations of 
as{Q), for Nf = 5 and Ap-^ chosen such that, in each order, q;s(Mzo) = 0.119. 

The functional form of as{Q), for 4 different values of Aj-j^ between 50 MeV and 350 MeV, is diplayed 
in figure m^a). The slope and dependence on the actual value of A;^ is especially pronounced at small 
Q^, while at large both the energy dependence and the dependence of A;^ becomes increasingly 
feeble. Nevertheless, as decreases with increasing energy scale and tends to zero at asymptotically high 
Q^. This is demonstrated in figure |H^b), where as(Q^) is plotted as a function of 1/ logQ. 

Any experimental proof of asymptotic freedom will therefore require precise determinations of ctg, 
or of other observables which depend on as(Q^), in a possibly large range of energy scales. This range 
should include as small as possible energies, since the relative energy dependence is largest there. To 
date, precise experimental data and respective QCD analyses are available in the range of Q ~ 1 GeV 
to a few hundred GeV. 

3.5 Relative size of finite order approximations 

The importance of higher order loop corrections and the degree of convergence of the perturbative 
expansion for the running ctg, see equation [71 is analysed and demonstrated in figure [HI In part (a), the 
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1-, the 2- and the 3-loop approximation of equation [71 each with A = 0.220 GeV, is shown. As can be 
seen, there is an almost 15% decrease of as when changing from 1-loop to 2-loop approximation, for 
the same value of A. The difference between the 2-loop and the 3-loop prediction is only about 1-2%, 
and it is less than 0.01% between the 3-loop and the 4-loop presentation which cannot be resolved in 
the figure. 

/ {4 — loop) {n — loop) \ 

The fractional difference in the energy dependence of Og, — alioov) y fo^' n = 1, 2 and 3, 

is presented in figure Efb). Here, in contrast to figure Efa), the values of A;^ were chosen such 
that Q!s(Mzo) = 0.119 in each order, i.e., A^g^ = 93 MeV (1-loop), A^ = 239 MeV (2-loop), and 
Aj^ = 220 MeV (3- and 4-loop). Only the 1-loop approximation shows sizeable differences of up to 
several per cent, in the energy and parameter range chosen, while the 2- and 3-loop approximation 
already reproduce the energy dependence of the best, i.e. 4-loop, prediction quite accurately. 

3.6 Quark masses and thresholds 

So far in this discussion, finite quark masses were neglected, assuming that both the physical and 
the renormalization scales and /x^, respectively, are larger than any other relevant energy or mass 
scale involved in the problem. This is, however, not entirely correct, since there are several QCD studies 
and «s determinations at energy scales around the charm- and bottom-quark masses of about 1.5 and 

4.7 GeV, respectively. 

Finite quark masses may have two major effects on actual QCD studies: Firstly, quark masses will 
alter the perturbative predictions of observables TZ. While phase space effects which are introduced by 
massive quarks can often be studied using hadronization models and Monte Carlo simulation techniques, 
explicit quark mass corrections in higher than leading perturbative order are available only for very few 
observables. 

Secondly, any quark-mass dependence of 71 will add another term /^^|^^^ to equation [TJ which 
leads to energy-dependent, running quark masses, ■mq{Q'^), in a similar way as the running coupling 
as(Q^) was obtained. 

In addition to these effects, as indirectly also depends on the quark masses, through the dependence 
of the (3 coefficients on the effective number of quarks flavours, Nf, with rUg <C /i. Constructing an 
effective theory for, say, (A^/-l) quark flavours which must be consistent with the Nf quark flavours 
theory at the heavy quark threshold fx^'^f^ ~ 0{mg), results in matching conditions for the ctg values of 
the (A^/-l)- and the A^/-quark flavours theories [SZj. 

In leading and in next-to- leading order, the matching condition is ai^^ = a^^ . In higher orders 
and the MS scheme, however, nontrivial matching conditions apply |3Z1 EHl ESj • Formally these are, if 
the energy evolution of as is performed in n^^ order, of order [n — 1). 

The matching scale fi^^f"^ can be chosen in terms of the (running) MS mass mg(/ig), or of the 
constant, so-called pole mass Mg. For both cases, the relevant matching conditions are given in |36j . 
These expressions have a particluarly simple form for the choice^ fx^^f"^ = mg{mq) or /i^^-f^ = Mg. In 
this report, the latter choice will be used to perform 3-loop matching at the heavy quark pole masses, 
in which case the matching condition reads, with a = ai^^'^ /n and a' = ai^^ /it: 

- = 1 + C2 + C3 a^ (8) 
a 

where C2 = -0.291667 and C3 = -5.32389 + {Nf - 1) ■ 0.26247 

The 4-loop prediction for the running as, using equation [7| with A^^"^"* = 220 MeV and 3-loop 

matching at the charm- and bottom-quark pole masses, ^i^^ = Mc = 1.5 GeV and = Mb = 

^The results of reference 1201 are also valid for other relations between fi'-^f'^ and rrig or Mg, as e.g. /i^^^^ = 2Mg. For 
3-loop matching, however, practical differences due to the freedom of this choice are negligible. 
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Figure 6: (a) 4-loop running of as{Q) with 3-loop quark threshold matching according to equations [3 
andlHl with A^^"^'' = 220 MeV and charm- and bottom-quark thresholds at the pole masses, /ii^-^"^'* = 

Mc = 1.5 GeV and fil^^~^^ = Mi, = 4.7 GeV (full line), compared with the unmatched 4-loop result 
(dashed line), (b) The fractional difference between the two curves in (a). 

4.7 GeV, is illustrated in figure (full line). Small discontinuities at the quark thresholds can be seen, 
such that ai^^ < ai^^^ by about 2 per mille at the bottom- and about 1 per cent at the charm- 
quark threshold. The corresponding values of Ajjg are A^^"^'' = 305 MeV and A^^"^'' = 346 MeV. 

Comparison with the 4-loop prediction, without applying threshold matching and for A^^^"^^ = 220 MeV 
and Nf = 5 throughout (dashed line) demonstrates that, in addition to the discontinuities, the matched 
calculation shows a steeper rise towards smaller energies because of the larger values of A^^"'^'* and 

a{^=3) 
MS 

The size of discontinuities and the changes of slopes are demonstrated in figure Eb? where the 
fractional difference between the two curves from figureEt, i.e. between the matched and the unmatched 
calculation, is presented. Note that the step function of ag is not an effect which can be measured; 
the steps are artifacts of the truncated perturbation theory and the requirement that predictions for 
observables at energy scales around the matching point must be consistent and independent of the two 
possible choices of (neighbouring) values of Nf. 

3.7 Perturbative predictions of physical quantities 

In practice, ctg is not an 'observable' by itself. Values of as(;U^) are determined from measurements of 
observables TZ for which QCD predictions exist. In perturbative QCD, these are usually given by a 
power series in as(/i^), like 
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= Pi(i?o + i?ias(;u2) + i?2(QV/^')«s(/^') + -) , 



(9) 



where Rn are the rith order coefficients of the perturbation series and PiRq denotes the lowest-order 
value of TZ. 

For processes which involve gluons already in lowest order perturbation theory, Pi itself may include 
(powers of) as- For instance, this happens in case of the hadronic decay width of heavy Quarkonia, 
r(T ggg hadrons) for which Pi cx a^. If no gluons are involved in lowest order, as e.g. in 
e"'"e~ ^ qq ^ hadrons or in deep inelastic scattering processes, PiRq is a constant and the usual 
choice of normalisation is P; = 1. Rq is called the lowest orc?er coefficient and Ri is the leading order 
(LO) coefficient. Following this naming convention, R2 is the next-to-leading order (NLO) and P3 is 
the next-to-next-to-leading order (NNLO) coefficient. 

QCD calculations in NLO perturbation theory are available for many observables TZ in high energy 
particle reactions like hadronic event shapes, jet production rates, scaling violations of structure func- 
tions. Calculations including the complete NNLO are available for some totally inclusive quantities, 
like the total hadronic cross section in e+e~ hadrons, moments and sum rules of structure functions 
in deep inelastic scattering processes, the hadronic decay widths of the Z° boson, of the r lepton and of 
heavy quarkonia like the T and the J /'^ . The complicated nature of QCD, due to the process of gluon 
self-coupling and the resulting large number of Feynman diagrams in higher orders of perturbation 
theory, so far limited the number of QCD calculations in complete NNLO. 

Another approach to calculating higher order corrections is based on the resummation of logarithms 
which arise from soft and collinear singularities in gluon emission j^lj. Application of resummation 
techniques and appropriate matching with fixed-order calculations are further detailed e.g. in 



3.8 Renormalization scale dependence. 

The principal independence of a physical observable TZ from the choice of the renormalization scale /i 
was expressed in equation [T] Replacing by as{^i^), using equation |21 and inserting the perturbative 
expansion of TZ (equation EI) into equation [T] results, for processes with constant P/, in 







+ 0{at) 



/i 



>dR2 



RiPo 



- [RiPi + 2P2^o 



(10) 



Solving this relation requires that the coefficients of a^ifi"^) vanish for each order n. With an appropriate 
choice of integration limits one thus obtains 



R, ^ 



Rq 


= const. 


Pi 


= const. 




= i?2(l) 








= ^3(1) 





Q2 

l3oRi In — , 



P3(l)-[2P2(l)/3o + Pi/3i]ln 



I + R,Pl In^ % 



as a solution of equation 1101 

Invariance of the complete perturbation series against the choice of the renormalization scale 
therefore implies that the coefficients Rn, except Pq and Pi, explicitly depend on /i^. In infinite order, 
the renormalization scale dependence of as and of the coefficients Rn cancel; in any finite (truncated) 
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order, however, the cancellation is not perfect, such that all realistic perturbative QCD predictions 
include an explicit dependence on the choice of the renormalization scale. 

This dependence is most pronounced in leading order QCD because Ri does not explicitly depend 
on fi and thus, there is no cancellation of the (logarithmic) scale dependence of as(/i^) at all. Only in 
next-to- leading and higher orders, the scale dependence of the coefficients Rn, for n >2, partly cancels 
that of asifi"^). In general, the degree of cancellation improves with the inclusion of higher orders in the 
perturbation series of TZ. 

Renormalization scale dependence is often used to test and specify uncertainties of theoretical calcu- 
lations and predictions of physical observables. In most studies, the central value of asifi"^) is determined 
or taken for /i equalling the typical energy of the underlying scattering reaction, like e.g. /i^ = E'^^ 
in e"'"e~ annihilation, and changes of the result when varying this definition of /i within "reasonable 
ranges" are taken as systematic uncertainties. 

There are several proposals of how to optimize or fix the renormalisation scale, see e.g. j^l^H^ I^ 
Unfortunately, there is no common agreement of how to optimize the choice of scales or how to define 
the size of the corresponding uncertainties. This unfortunate situation should be kept in mind when 
comparing and summarising results from different analyses. 

For more details and examples on renormalisation scale dependences see e.g. 

3.9 Nonperturbative QCD methods 

At large distances or low momentum transfers, as becomes large and application of perturbation theory 
becomes inappropriate. Nonperturbative methods have therefore been developed to describe strong 
interaction processes at low energy scales of typically < 1 GeV^, such as the fragmentation of quarks 
and gluons into hadrons ( "hadronisation" ) and the absolute masses and mass splittings of mesons. 

Hadronisation models are used in Monte Carlo approaches to describe the transition of quarks and 
gluons into hadrons. They are based on QCD-inspired mechanisms like the "string fragmentation" 
ES] or "cluster fragmentation" [IHl- Those models usually contain a number of free parameters 
which must be adjusted in order to reproduce experimental data well. They are indispensable tools not 
only for detailed QCD studies of high energy collision reactions, but are also important to assess the 
resolution and acceptance of experimental setups. 

Power corrections are an analytic approach to approximate nonperturbative hadronisation effects 
by means of perturbative methods, introducing a universal, non-perturbative parameter 



to parametrize the unknown behaviour of as{Q) below a certain infrared matching scale /x/ |17j. Power 
corrections are regarded as an alternative approach to describe hadronisation effects on event shape 
distributions, instead of using phenomenological hadronisation models. 

Lattice Gauge Theory is one of the most developed nonperturbative methods (see e.g. jlSlESl) and 
is used to calculate, for instance, hadron masses, mass splittings and QCD matrix elements. In Lattice 
QCD, field operators are applied on a discrete, 4-dimensional Euclidean space-time of hypercubes with 
side length a. 

4 Tests of Asymptotic freedom in e+e~ annihilations 

While historical and early developments of QCD were mainly inspired by the results from deep inelastic 
lepton-nucleon scattering experiments, hadronic final state of e"'"e~ annihilations developed into a prime- 
tool for precision tests of QCD, due to the availability of higher effective energies and the point-like 
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Figure 7: Energy dependence of hadronisation corrections to 3-jet event production rates, defined using 
various jet reconstruction schemes. 

and simple nature of the primary particle reaction. Studies of hadron jets and determinations of as, 
in a wide range of centre of mass energies and with increasing experimental and theoretical precision, 
led to many signatures of asymptotic freedom and, equivalently important, to the evidence for colour 
charged gluons and the SU(3) gauge structure of QCD. 

4.1 Jet production rates 

The study of multijet event production rates in e+e~ hadronic final states gave direct evidence for 
the energy dependence of ag and for asymptotic freedom, before summaries of explicite measurements 
of as were able to provide a convincing case (c.f. figure Ej). The basic idea was rather simple and 
straight-forward. It was based on a definition of resolvable jets of hadrons which could be applied both 
in perturbative QCD calculations and on measured hadronic final states: 

Within the JADE jet algorithm |j25j, the scaled pair mass of two resolvable jets i and j, yij = 
Mfj/E^j^g, is required to exceed a threshold value ycut, where Eyis is the sum of the measured energies 
of all particles of a hadronic final state - or, in theoretical calculations, of all quarks and gluons. In a 
recursive process, the pair of particles or clusters of particles n and m with the smallest value of t/nm 
is replaced by (or "recombined" into) a single jet or cluster k with four-momentum Pk = Pn + Pm, as 
long as ynm < Vent- The procedure is repeated until all pair masses yij are larger than the jet resolution 
parameter T/cut, and the remaining clusters of particles are called jets. 

Several jet recombination schemes and definitions of Mij exist [M.; the original JADE scheme with 
Mfj = 2 ■ Ei ■ Ej ■ {1 — cos 9ij), where Ei and Ej are the energies of the particles and 9ij is the angle 
between them, and the "Durham" scheme with Mfj = 2 • mm{Ef,Ej) ■ (1 — cos%), were most 
widely used at LEP, due to their superior features like small sensitivity to hadronisation and particle 
mass effects jlHl- 

The virtues of the JADE jet definition were especially suited for an early observation of the energy 
dependence of a^, without the need to determine ag and therefore avoiding the large systematic and - 
at that time - theoretical uncertainties: 

• The relative production rate R3 of 3-jet hadronic final states follows a particular simple theoretical 
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expression: 

^3 = ^f^, ^ J^^^l = Ci (?/cut) as(Ai^) + C2(?/cut,/i^)as(/i^) , 

cj(e+e~ nadrons) 

in NLO perturbation theory. In leading order, is thus directly proportional to ag, and any 
energy dependence of R3 observed in the data must be due to the energy dependence of - if there 
are no other energy dependent effects. The coefficients Ci and C2 are energy independent. They 
can be reliably calculated and predicted by QCD, whereby the renormalisation scale dependence 
of C2 is only a small disturbance. 

Model studies showed that hadronisation corrections to are small and, in a suitable range of 
centre of mass energies, almost constant, see figure [7| pHj . 

The JADE jet algorithm is particlarly easy to apply to measured hadronic final states, and cor- 
rections due to limited detector resolution and acceptance are small and manageable. 
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Figure 8: Energy dependence of 3-jet event production rates, measured using the JADE jet finder at 
a scaled jet energy resolution y^ut = 0.008. The errors are experimental. The data are not corrected 
for hadronisation effects. They are compared to theoretical expectations of QCD, of an abelian vector 
gluon model, and to the hypothesis of a constant coupling strength. 

The first experimental study of the energy dependence of 3-jet event production rates, at cm. 
energies etween 22 and 46 GeV, analysed for constant jet resolution ?/cut at the e+e~ collider PETRA, 
gave first evidence for the energy dependence of as already in 1988 [22] • These data are shown in 
figure IHl together with more results from eperiments at the PEP, TRISTAN [21] and finally, at the 
LEP collider |S2]- The measured 3-jet rates significantly decrease with increasing centre of mass energy, 
in excellent agreement with the decrease predicted by QCD. The hypthesis of an energy independent 
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coupling, and especially the prediction of an alternative, QED-like abelian vector gluon model, where 
gluons carry no colour charge, are in apparent contradiction with the data 

In order to further demonstrate asymptotic freedom with these data, they are - combined at suitable 
mean energies - plotted against 1/ In Ecm, as shown in figure IHl For infinite energies, Ecm oo, as and 
thus i?3 are expected to vanish to zero, which is in very good agreement with the data. 



E,„[GeV] 



© © ^ 

© © OS 



CO 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.05 0.10 0.15 0.20 0.25 0.30 0.35 

1 / ln(E,„/GeV) 



Figure 9: 3-jet event production rates as shown in Fig. |H1 however as a function of l/lnEcm, to 
demonstrate that i?3 oc ctg — ^ at asymptotic (i.e. infinite) energies. 



4.2 Evidence for the gluon self coupling 

The gluon self-coupling, as a direct consequence of gluons carrying colour charge by themselves, is 
essential for the prediction of asymptotic freedom. A rather direct method to detect effects of gluon- 
selfcoupling was accomplished at the LEP collider, by analysing distributions which are sensitive to 
the spin structure of hadronic 4-jet final states j2ZI- For instance, the so-called Bengtson-Zerwas angle, 
Xbz [53], measuring the angle between the planes defined by the two highest and the two lowest energy 
jets, is rather sensitive to the difference of a gluon-jet splitting into two gluons, which in QCD is the 
dominant source of 4-jet final states, and a gluon splitting to a quark-antiquark pair, which is the 
dominant process in an abelian vector theory where gluons carry no colour charge. 

The results of an early study which showed convincing evidence for the gluon self coupling [S^ after 
only one year of data taking, is shown in figure ^| The data clearly favour the QCD prediction and 
rule out the abelian vector gluon case. 
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Figure 10: Distribution of the azimuthal angle between two planes spanned by the two high- and the 
two low-energy jets of hadronic 4-jet events measured at LEP |54j, compared to the predictions of QCD 
and of an abelian vector gluon model where gluons carry no colour charge |27j . 

4.3 Determination of the QCD group constants 

As already mentioned in section 3, the QCD group constants Ca, Cp and Nj, assume values of 3, 4/3 
and 5 in a theory exhibiting SU(3) symmetry, with 5 quark flavours in the vacuum polarisation loops. 
For alternative, QED-like toy models of the strong interaction with U(l) symmetry, these values would 
be Ca = and Cp = 0.5. Experimental determination of these two parameters can thus be regarded 
as one of the most intimate tests of the predictions of QCD. 

At LEP, data statistics and precision allowed to actually determine experimental values for Ca, 
which basically is the number of colour charges, and Cp. The current state-of-the art of such studies, 
involving analyses of several 4-jet angular correlations or fits to hadronic event shapes, is summarised 
j55j in figure ^3 The data, with combined values of 

Ca = 2.89 ± 0.01 (stat.) ± 0.21 (syst.) (12) 
Cp = 1.30 ± 0.01 (stat.) ± 0.09 (syst.) 

are in excellent agreement with the gauge structure constants of QCD. They rule out the Abelian vector 
gluon model and theories exhibiting symmetries other than SU(3). 

4.4 The running as from 4-jet event production 

Since the beginning of the 1990's, advances in phenomenological calculations and predictions as well as in 
experimental techniques allowed to perform precision determinations of as, in a broad range of energies 
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Figure 11: Measurements and combination of the QCD colour factors Ca and Cp 



and based on many different methods and observables. While an overall summary of measurements, 
as presented in section 6, gives a very distinct signature for asymptotic freedom, the demonstration 
of the running of from single experiments, minimising point-to-point systematic uncertainties, adds 
extra confidence in the overall conclusion. 

One of the most recent developments in determining in e+e^ annihilation is the precise extraction 
of from differential 4-jet distributions, i.e. the 4-jet production rate as a function of the jet resolution 
Vent inni inn IHH] • QCD predictions in NLO are available [221, which is in 0[al) perturbation theory. 
4-jet final states appear only at 0{a1) as their leading order, since they require at least two radiation 
or splitting processes off the primary quark- ant iquark pair. Although "only" in NLO, the theoretical 
uncertainties for 4-jet observables appear to be rather small, smaller than for typical 3-jet event measures 
or event shapes. Since 4-jet observables are proportional to , they provide very sensitive measurements 
of a^. 

As another new development, the data of the JADE experiment, at the previous PETRA collider 
which was shut down in 1986, are currently re-analysed, with the experimental and theoretical experience 
and tools of today. This adds important new information especially at lower energies, seen from the 
perspective of the LEP experiments. Moreover, the JADE experiment can be regarded as the smaller 
brother or prototype of the OPAL detector at LEP, since both experiments were based on similar 
detector techniques. 

A LEP-like study of from 4-jet production from JADE is available jSOl, and is summarised in 
figure IT^ together with the corresponding results from OPAL and the other LEP experiments. In order 
to appreciate the high significance for the running of cts, one should keep in mind that only the inner 
(experimental) error bars need to be considered in a relative comparison of these data. 
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Figure 12: The values of ag from 4-jet event producion. Errors are experimental (inner marks) and 
the total errors IHOl EZl EH]- The lines indicate the QCD prediction for the running of as with 
a,(Mzo) = 0.1182 ± 0.0027 

5 QCD tests in deep inelastic lepton-nucleon scattering 

As outlined in section 2 and 3, the observation of approximate scaling of nuclear structure functions, and 
thereafter, with higher precision and extended ranges of x and Q^, of (logarithmic) scaling violations, 
originally boosted the development of the quark-parton model and of QCD. The limited range of 
fixed-target lepton-nucleon scattering experiments in x and Q^, however, prevented significant and 
unambiguous tests of QCD scaling violations and the running of ctg, see e.g. 

This picture changed dramatically when the HERA electron-proton and positron-proton collider 
started operation in 1991, with lepton beam energies of 30 GeV and protons of 920 GeV. HERA 
extended the range in by more than 2 orders of magnitude towards higher values, and the range 
in X by more than 3 orders of magnitude towards smaller values. With these parameters, precise tests 
of scaling violations of structure functions, but also precise determinations of the running as from 
jet production were achieved. While these two topics will be reviewd in the following subsections, a 
summary of significant determinion in deep inelastic scattering will be included in section 6, see also 

El EH EH]- 

5.1 Basic introduction to structure functions 

Cross sections of physical processes in lepton-nucleon scattering and in hadron-hadron collisions depend 
on the quark- and gluon-densities in the nucleon. Assuming factorisation between short- distance, hard 
scattering processes which can be calculated using QCD perturbation theory, and low-energy or long- 
range effects which are not accessible by perturbative methods, such cross sections are parametrized 
by a set of structure functions Fi {i= 1,2,3). The transition between the long- and the short-range 
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regimes is defined by an arbitrary factorization scale /x/, wliicli — in general — is independent from 
tlie renormalization scale n, but has similar features as the latter: the higher order coefficients of the 
perturbative QCD series for physical cross sections depend on /ij in such a way that the cross section 
to all orders must be independent of /x/, i.e. ^jda/dfif = 0. To simplify application of theory to 
experimental measurements, the assumption fif = fi is usually made, with /i = Q as the standard 
choice of scales. 

In the naive quark-parton model, i.e. neglecting gluons and QCD, and for zero proton mass, the 
differential cross sections for electromagnetic charged lepton (electron or muon) - proton scattering off 
an unpolarized proton target is written 



dV- Arra' L ,2,^em , - V) , ^em 



[1 + (1 + + ^ -{F^"" - 2xF^ 



X 



(13) 



where, for fixed target reactions, x = 2M(e-e') nucleon momentum fraction carried by the struck 

parton, y = 1 — E' /E, is the negative quadratic momentum transfer in the scattering process, 
and M, E and E' are the mass of the proton and the lepton energies before and after the scattering, 
respectively, in the rest frame of the proton. In the quark-parton model, these structure functions 
consist of combinations of the quark- and antiquark densities q{x) and q{x) for both valence- (u,d) and 
sea-quarks (s,c). 

In QCD, the gluon content of the proton as well as higher order diagrams describing photon/gluon 
scattering, 7(7 qq, QCD Compton processes 7g — > gq and gluon radiation off quarks must be taken 
into account. Quark- and gluon-densities, the structure functions Fi and physical cross sections become 
energy (Q^) dependent. QCD thus predicts, departing from the naive quark-parton model, scaling 
violations in physical cross sections, which are associated with the radiation of gluons. While perturba- 
tive QCD cannot predict the functional form of parton densities and structure functions, their energy 
evolution is described by the so-called DGLAP equations [HT] IH^. 

The energy dependence of structure functions are known, since the early 1980's, in NLO QCD |63j . 
Recently, parts of the NNLO predictions became available, see e.g. jHSlElI, and should be completely 
known in NNLO fairly soon. Apart from terms whose energy dependence is given by perturbative QCD, 
structure functions contain so-called "higher twist" contributions (HT). The leading higher twist terms 
are proportional to l/Q^; they are numerically important at low < 0(few GeV^) and at very large 
a; ~ 1. 



5.2 Scaling violations of structure functions 

A recent summary of measurements of the proton structure function F2 is given in figure El j66, . The 
measurements cover an impressively large parameter space in x and Q^. The degree of scaling violations 
predicted by QCD ^T], namely a strong increase of with increasing at small x, and a decrease of 
F2 at X > 0.1, is clearly reproduced by the data. Over the whole kinematic range, the data are in very 
good agreement with the structure function evolution as predicted by QCD. 

One of the most important results obtained from HERA was, in fact, the discovery of the strong 
rise of F2 at small x, as demonstrated again in figure El It proves the QCD picture of an increasing 
part of the proton consisting of gluons, i.e. field energy, if the resolution is increased (c.f. section 2.5). 
Another way to demonstrate the compatibility of the observed scaling violations with the predictions of 
QCD is to analyse the lorgarithmic slopes, dF2{x, Q^)/dlnQ'^ as a function of x, as shown in figure IT31 
Again, QCD provides an excellent description of the data. 
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Figure 13: Summary of measurements of F2 [B^. For better visibility, the results for different values of 
X were multiplied with the given factors of 2*. 



5.3 The running as from jet production in DIS 

The increase of momentum transfer at HERA also allowed to see and analyse the jet structure in 
deep inelasting scattering processes. Here, similar as in e"'"e~ annihilations, production of multijet final 
states is predicted in QCD to NLO, and allows to determine as in a large range of energy scales. 

Inclusive as well as differential jet production rates were studied in the energy range of ~ 10 
up to 10000 GeV^, using similar jet definitions and algorithms as in e+e~ annihilation. In leading 
order Og, 2 + 1 jet events in deep inelastic ep scattering arise from photon-gluon fusion and from QCD 
Compton processes. The term '2 + 1 jet' denotes events where two resolved jets can be identified at 
large momentum transfer, in addition to the beam jet from the remnant of the incoming proton. 

A recent summary of as determinations from the two HERA experiments HI and ZEUS is given 
in figure El [HZ]. Here, the transverse jet energy Ej^* was chosen as the relevant energy scale. Both 
experiments have determined as at several different values of E^^^, and the summary of all these results 
clearly demonstrates that as runs as predicted by QCD. 

A combination of these results jHZI will be included in the overall summary of as determinations, 
which is presented in the following section. 
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Figure 14: Proton structure function F2{x,Q'^) [HE], as a function of x, demonstrating the strong rise 
at very small values of x. 

6 Summary of measurements 

6.1 Previous developments and current status 

By the end of the 1980's, as concluded in section 2.11 above, summaries of measurements of demon- 
strated very good agreement with the prediction of asymptotic freedom (c.f. figure E)), but both the 
precision and the energy range of the results did not yet allow to claim firm evidence for the running 

of ttg. 

This situation changed significantly with the turn-on of the LEP e"^e~ collider in 1989, with ad- 
vancements in theoretical calculations and predictions, and with improved experimental methods to 
determine as- At LEP, after the first successful year of data taking, precise results emerged at high 
energies, i.e. at the Z° mass peak, Ecm = Mzo = 91.2 GeV, and also at very small energy scales of 
Mr = 1.78 GeV, extending the energy range and improving the overall precision by significant amounts. 

This in turn quickly led to a very consistent picture: the experimental results agreed with the 
prediction of asymptotic freedom in the energy range from the r to the Z° mass, and were compatible 
with a common value of as(Mzo) = 0.117 ± 0.004 jUH]- The theoretical uncertainties of this average 
were, and still are today, a matter of discussion. 

The number of precise measurements of continued to grow, especially from the HERA collider, 
from hadron collisions at the Tevatron, and from the study of heavy quarkonia masses and decays. 
Through LEP-II and the Tevatron, the energy scale of determinations was extended to 200 GeV 
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Figure 15: Logarithmic slope of F2{x, Q"^) as a function of x . 

and beyond. The precision of many measurements and theoretical calculations was increased, and the 
number of observables and measurements which were based on complete NNLO QCD predictions grew. 
A comprehensive summary of the state-of-the-art of determinations was given e.g. in \^2^ and was 
continuously updated in jHHinni, leading to a combined world average of a^^Mio) = 0.1182 ± 0.0027 in 
20046. 

In the following, these previous summaries will be updated, the most recent and new results will be 
presented, the actual evidence for asymptotic freedom will be demonstrated and a new world average 
of as(Mzo) will be derived. 

6.2 Update: new results since 2004 

Since the last summary in 2004 i69j, a number of new results and updates of previous studies were 
presented. The most significant which will be included in the overall summary and determination of 
the world average of as(Mzo) were: 

• A re-evaluation jTUj of the existing data on hadronic decays of the r-lepton and a revision of 
the theoretical framework, especially in terms of narrowing the range of theoretical uncertainties, 
resulted in a significantly improved value of as(M^) = 0.345 ± 0.010, in full NNLO QCD. The 
total error is thereby reduced by a factor of 3, w.r.t. previous estimates [7T] - which is basically 
due to the inclusion of available terms in NNNLO and a rather restrictive treatment of systematic 

^For alternative summaries of ag, leading to almost identical values of as{Mzo) with a different treatment of the overall 
uncertainties, see e.g. j31j 



26 



HERA 



0.25 




T 



ZEUS (a) 
HI 



0.2 



0.15 




0.1 



QCD 

a/M^.) = 0.118 + 0.003 



10 




(GeV) 



Figure 16: Results of as as a function of Ej? from HERA experiments HI and ZEUS jUTj . 

uncertainties. Although the error on this result, with a wider treatment of systematics, can well 
increase by a factor of 2 ^72j, the published value is retained for further analysis in this review. 

• New analyses and a new combination of results from jet production in deep inelastic electron or 
positron - proton scattering at HERA [HT], as shown in figure 16, provided an improved overall 
value of as(Mzo) = 0.1186 ± 0.0051, in NLO of perturbative QCD. 

• A new study of hadron masses using predictions from lattice gauge theory, including vacuum polar- 
isation effects from all three light quark flavours and improved third and higher order perturbative 
terms, resulted in a new and improved value of as(Mzo) = 0.1170 ± 0.0012 [73]. Although the 
methods used in this study and the small size of the claimed overall error are still under discussion 
[7Ij . the published value is retained here for further discussion. 

• New studies of 4-jet final states in e"'"e~ annihilation at LEP |56t 15 7| 13^ I60j . see also section 4.4, 
and a combination of the respective as results give a new average of as(Mzo) = 0.1176 ± 0.0022, 
in 0{al) which, for 4-jet production, corresponds to NLO in perturbative QCD. 

In the following overall summary of measurements of a^, these four results will replace the respective 
values used in the previous summary of 2004 [UUj . 

6.3 as summary 

The new overall summary of as is given in table ^ where the new and updated results discussed in 
the previous section are underlined. Most of the results given in table Q are combined from several 
or many individual measurements of different experiments and groups. For results obtained at fixed 
energy scales Q (or in narrow ranges of Q), the value of as{Q) is given, together with the extrapolation 
to the "standard" energy scale, Q = Mgo, using equation [7| in 4-loop approximation and 3- loop quark 
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threshold matching at the heavy quark pole masses Mc = 1.5 GeV and Mb = 4.7 GeV. Results from 
data in ranges of energies are only given for Q = Mzo. Where available, the table also contains the 
contributions of experimental and theoretical uncertainties to the total errors in as(Mzo). 

Finally, in the last two columns of table the underlying theoretical calculation for each mea- 
surement and a reference to this result are given, where NLO stands for next-to-leading order, NNLO 
for next-next-to-leading-order of perturbation theory, "resum" stands for resummend NLO calculations 
which include NLO plus resummation of all leading und next-to-leading logarithms to all orders (see 
and and "LGT" indicates lattice gauge theory. 
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Figure 17: . Summary of measurements of as{Q) as a function of the respective energy scale Q, from 
tabled Open symbols indicate (resummed) NLO, and filled symbols NNLO QCD calculations used in 
the respective analysis. The curves are the QCD predictions for the combined world average value of 
as(Mzo), in 4-loop approximation and using 3-loop threshold matching at the heavy quark pole masses 
Mc = 1.5 GeV and = 4.7 GeV. 

In figure [T7[ all results of as{Q) given in table [T] are graphically displayed, as a function of the 
energy scale Q. Those results obtained in ranges of Q and given, in table [T[ as a^iMzo) only, are not 
included in this figure - with one exception: the results from jet production in deep inelastic scattering 
are represented in tabled by one line, averaging over a range in Q from 6 to 100 GeV, while in figure IT7I 
combined results for fixed values of Q as presented in are displayed. 
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The data are compared with the QCD prediction for the running as, calculated for as(Mzo) = 
0.1189 ±0.0010 which - as will be discussed in the next subsection - is the new world average. The QCD 
curves are calculated using the 4-loop perturbative prediction, equation and 3-loop quark threshold 
matching, see section 3.6. The data are in excellent agreement with the QCD prediction, from the 
smallest to the largest energy scales probed by experimental data. 
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Figure 18: Some of the most significant data from figure El however plotted as a function of 
1/ log(Q/GeV^), in order to demonstrate the significance of the data showing evidence for asymptotic 
freedom, i.e. the vanishing of as{Q) at asymptotically large energy scales. 

In order to demonstrate the significance with which the data probe asymptotic freedom, some of the 
most precise ag results are plotted^ in figure now as a function of the inverse logarithm of Q. Also 
shown in this figure is the QCD prediction in 4-loop approximation including 3-loop quark threshold 

^Only a selection of results, especially of those at high energy scales, is shown in this figure, to maintain visibility and 
clarity. The measurements chosen for this demonstration, from left to right, are: as from e~''e~ event shapes and jets at 
LEP-II, for {Q =189 GeV), from 4-jet events at LEP (Q = Mz" — 91.2 GeV), from the reanalysis of e"''e~ event shapes 
and jets at PETRA {Q = 44 GeV), from hadron masses and lattice theory a,t Q — 7.5 GeV, from heavy quarkonia decays 
at Q = 4.75 GeV, and for r-decays at Q = 1.78 GeV. 
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matching, for as{Mzo) = 0.1189 (full line), and, for demonstration only, the 4-loop QCD curve omitting 
quark threshold matching. 

The data very significantly prove the particular QCD prediction of asymptotic freedom. Apart from 
precisely reproducing the characteristic QCD-shape with an inverse logarithmic slope, the data point 
at the very lowest energies, i.e. the right-most point in figure ^1 indicates that the available precision 
allows to conclude that quark threshold matching is necessary for QCD to consistently describe the 
data. 

In fact, data precision is now so advanced that a rather simple QCD fit, e.g. in leading order QCD 
with no threshold matching, with a fit probability of less than 1%, fails to describe data^. Evidently, the 
probabilities for a hypothetical constant and energy mdependent ag^, or an abelian vector gluon theory 
which predicts an increase of the coupling with increasing energy scale, c.f. figure |Hl have negligible 
probabilities to describe data. The same is true for other functional forms, like oc 1/Q or c^s oc 1/Q^ 
- such functions may be adjusted such that they can describe a few data points either at low or at high 
values of Q, but altogether fail to describe data in the full range of energy scales from 1.78 to 200 GeV. 

Therefore, it is concluded that the data, with the current precision which has substantially increased 
over the past few years, prove the specific QCD functional form of the running coupling as, and therefore 
of asymptotic freedom. 

6.4 A new world average of as{Mzo) 

Since all measurements of as, as shown in figure IT71 are consistent with the running coupling as predicted 
by QCD, it is legitimate to evolve all values of as{Q) to a common reference energy scale, like e.g. 
Q = Mio, using the QCD 4-loop beta-function. This was already done in the past, see e.g. [321 1^ - 
and world average values of as(Mzo) were determined. 

A general difficulty in the averaging procedure is the proper and optimal treatment of errors and 
uncertainties. For many of the results, the dominating error is the theoretical uncertainty - which 
however is not uniquely defined, and is treated quite differently in most of the studies. While in a 
majority of analyses the theoretical uncertainty is defined through a variation of the renormalisation 
scale, the range of this variation is not unique. Theoretical errors, in some studies, are defined and 
probed in a rather restrictive and optimistic way, while in other studies they are treated in a more 
generous manner. Also, theoretical uncertainties are likely to be correlated between different analyses, 
however to an unknown degree. So the assigned errors cannot be treated like statistical and independent 
errors. 

In previous studies [HSlEIllSni, an error- weighted average and an "optimised correlation" was cal- 
culated from the error covariance matrix, assuming an overall correlation factor between the errors of 
the measurements included in the averaging process. This factor was adjusted such that the overall 
equals unity per degree of freedom [7S]. Usually the overall was larger than unity per d.o.f., which 
indicated that either the individual errors were overestimated, or they were correlated. Adjusting an 
overall correlation factor should take care about the latter case. 

The procedure, applied in a previous study to all results which are based on NNLO QCD and which 
have overall errors of < 0.008 on as(Mzo), then resulted in as(Mzo) = 0.1182 ± 0.0027 Here, the 
averaging shall be done using the same method, however the selection of results to be included in the 
final average shall be altered: 

First, the restriction to results which are complete to NNLO will give a strong weight to studies at 
low energy scales, because most of the results obtained at high energies are in NLO, possibly including 
resummation, only. Second, it may be time to include the results from lattice gauge calculations, which 

^Such a "simple" fit was previously used to "fit" the /?o coefficient of the QCD beta-function, c.f. equation|31 or - 
alternatively - the number of colour degrees of freedom, Ca — = i- 
^In fact, there exists no theory which predicts a constant coupling. 
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now have reached a maturity which may be comparable to those in (NNLO) perturbation theory - and 
lattice gauge results are most likely not correlated with classical perturbation theory, thereby adding an 
independent input to the standard selection. Finally, in order to reach a higher degree of independence 
between the results included in the averaging procedure, data from different processes should be chosen, 
avoiding a possible bias towards any direction. In this sense, the precise determinations of as from 4-jet 
events in e^e~ annihilation at Q = Mz, from event shapes and jets at Q = 189 GeV and from jets in 
deep inelastic scatterin are also chosen to be included. 

Follwing this strategy, 10 results are henceforth selected to be included in the averaging procedure. 
They are summarised in table|21and average to as{Mzo) = 0.1189 ±0.0007 , with an overall of 9.9 for 
9 degrees of freedom. Since the is larger than unity per d.o.f., no common correlation factor needs 
to be assumed; in fact, in order to reach exactly 1 per d.o.f., the assigned errors of single measurements 
should be increased, a method which is frequently used e.g. in [21]. 

The fact that - in contrast to the case of previous reviews, see e.g. [^211221 - I d.o.f. is not smaller 
than unity is mainly caused by two of the updates which were discussed above, namely by the new 
assessment of r decays and from heavy hadron masses in lattice theory. Both have the smallest overall 
errors assigned - ±0.0012 on as{Mio). Treating these as gaussian and independent errors, the two 
results are 2.7 standard deviations apart from each other. 

The question whether one or both of these measurements has underestimated its assigned overall 
uncertainty cannot finally be answered. The significance of each of the 10 selected measurements can 
be judged from the last column of table 2, which shows the deviation of the respective average value 
of ag(Mzo) when omitting this particular measurement in the averaging procedure. The maximum 
deviation observed is ±0.0011 — 0.0013, a value which is compatible with a total error of ±0.0010. 
The largest change in is observed when the result from r decays is left out - a possible hint for an 
underestimation of its assigned error. 

Leaving out two of the 10 selected measurements when averaging the results gives values of as(Mzo) 
which vary between 0.1173 and 0.1205; these two extremes again average to 0.1189 with a maximum 
deviation of 0.0016. 

In view of these studies and variations, it is finally concluded that 



a,(Mzo) = 0.1189 ±0.0010 



is the new world average^° of a^. Here, the overall error decreased by almost a factor of three as compared 
to the previous review jBHI- This small error, however, appears to be realistic since all measurements 
agree well with this new average, as can be seen in figure IT7I the error band is very narrow, but all the 
data are consistent with this result. 

We have therefore reached, after 30 years of QCD, the case that not only asymptotic freedom is 
proven, beyond any doubts, by the data, but also that as{Mio) is now known very precisely, to better 
than 1% accuracy! 

7 Summary and Outlook 

The concept of asymptotic freedom, i.e. the QCD prediction of an inverse logarithmic decrease of the 
coupling strength with the energy or the momentum transfer in high energy scattering reactions, was 
shown to be significantly and reliably verified by a number of different measurements. Historically, the 
first signatures for asymptotic freedom came from the observation of approximate scaling of structure 
functions in deep inelastic lepton-nucleon scattering experiments, and the subsequent observation of 

-"^"^A small increase of the error of 0.0007 from the covariance matrix seems justified, due to the scatter of averages when 
leaving out one or two of the results, and due to being slightly larger than 1 per d.o.f. 
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small scaling violations, in a manner as was predicted by QCD. Before specific measurements of as at 
different energy scales and processes were precise enough to proove the running of as, measurements 
of the energy dependence of jet production rates in e"'"e~ annihilations provided strong evidence for 
the energy dependence of as, as predicted by QCD and by the concept of asymptotic freedom. The 
process of gluon self-coupling was experimentally established by studies of angular correlations within 
4-jet final states, and the SU(3) gauge structure of QCD was confirmed by studies of 4-jet and event 
shape observables at the LEP e^e~ collider. 

Measurements of as, performed by single experiments over ranges of energy scales, demonstrate 
evidence for the running of as, in perfect agreement with QCD. Studies by single (or by few but similar) 
experiments avoid systematic uncertainties which vary from (energy) point to point, and therefore have 
the potential to demonstrate asymptotic freedom, however, in most cases only over a limited range of 
energies. Examples of such studies, as as from jet production in deep inelastic lepton-nucleon scattering, 
as from 4-jet event production in e^e~ annihilation, are presented and discussed in this review. 

The most significant experimental proof of asymptotic freedom today is provided by the summary 
and combination of all measurements of as, over an energy range of 1.6 GeV to more than 200 GeV, 
from all available processes and experiments, involving perturbative and lattice QCD calculations. The 
results are in excellent agreement with QCD and precisely reproduce the inverse logarithmic dependence 
of as from the energy or momentum transfer scale Q. The data proove the neccessity to include higher 
loop diagrams in the perturbative prediction of the energy dependence of as, and demonstrate the need 
for quark flavour threshold matching of as, again as predicted by QCD. 

Since all measurements of as are in excellent agreement with asymptotic freedom as predicted by 
QCD, the results are elvolved, according to the QCD 4-loop beta-function, to a common energy scale 
Q = Mip. A set of significant and precise measurements, well balanced over all available processe, 
energies, experiments and theoretical methods, result in a new and improved world average value of 

a^(Mzo) = 0.1189 ±0.0010 . 

The overall uncertainty of this result is improved by almost a factor of 3, compared to the previous 
average presented in 2004. 

This improvement will have significant implications on verifying the grand unification of forces, and 
in particular, on possible signatures for Supersymmetry see e.g. jUUIHSl- The improvement is achieved 
by the inclusion of several new results and updates, like the one from r decays, from 4-jet final states 
in e"'"e~ annihilations at LEP, and from lattice predictions of heavy hadron masses. 

The total uncertainties of these new results are in the range of 1 to 2 % only, and therefore these 
measurements largely dominate the determination of the world average. Although some of the quoted 
uncertainties may be underestimated, no significant disagreement between any of the measurements nor 
with the world average value of as(Mzo) was found. 

It is therefore concluded that the running coupling averages to the value of as(Mzo) given above, 
with an overall uncertainty now being less than 1%. 

Future prospects for precise determinations of as are upcoming QCD calculations, in NNLO pertur- 
bation theory, for more observables and processes, which will replace the currently used calculations in 
NLO or in resummed NLO, see e.g. tabled Such calculations are in preparation by several groups, see 
e.g. inni inni ISl] , and the results are eagerly awaited by experimentalists. It is hoped that application 
of these new NNLO calculations will provide an increased number of precise as determinations, with 
overall errors of the order of 1%. This will provide the means for precise compatibility checks and for 
further improvements of the world average value of as. 
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Table 1: World summary of measurements of as (status of April 2006): DIS = deep inelastic scattering; 
GLS-SR = Gross-Llewellyn-Smith sum rule; Bj-SR = Bjorken sum rule; (N)NLO = (next-to-) next-to- 
leading order perturbation theory; LGT = lattice gauge theory; resum. = resummed NLO. New or 
updated entries since the review of 2004 [^2] are underlined. 





Q 






Aas(Mzo) 






Process 


[GeV] 


«.(g) 


«s(Mzo) 


exp. 


theor. 


Theory 


refs. 


DIS [pol. SF] 


0.7 - 8 




Q 1 1 + 0.010 
U.11-3 _ 0.008 


±0.004 


+0.009 
-0.006 


NLO 


m 


DIS [Bj-SR] 


1.58 


n Q7e; + 0.062 
U.0(0 _ 0,081 


n 1 01 + 0.005 

U.iZl „ 0.009 


- 


- 


NNLO 


m 


DIS [GLS-SR] 


1.73 


0.280 t 


n 1 1 + 0.009 

U.ilZ _ 0.012 


+0.008 
-0.010 


0.005 


NNLO 


m 


r-decays 


1.78 


0.345 ±0.010 


0.1215 ±0.0012 


0.0004 


0.0011 


NNLO 


m 


DIS [v- xFa] 


2.8 - 11 




^ . ^ Q + 0.007 

U-iiy - 0.006 


0.005 


+0.005 
-0.003 


NNLO 


m 


DIS [e//x; F2] 


2 - 15 




0.1166 ±0.0022 


0.0009 


0.0020 


NNLO 


IHOIIHI] 


DIS [e-p — i> jets] 


6 - 100 




0.1186 ±0.0051 


0.0011 


0.0050 


NLO 




T decays 


4.75 


0.217 ±0.021 


0.118 ±0.006 


- 


- 


NNLO 


m 


QQ states 


7.5 


0.1886 ±0.0032 


0.1170 ±0.0012 


0.0000 


0.0012 


LGT 


m 


e+e- [Fl] 


1.4 - 28 




n 1 1 QQ + 0.0044 

u.iiyg _ 0.0054 


0.0028 


+ 0.0034 
- 0.0046 


NLO 


m 


e+e- [cThad] 


10.52 


0.20 ±0.06 


1 30 + °-°2i 

U.iOU _ 0.029 


+ 0.021 
- 0.029 


0.002 


NNLO 


m 


e+e- [jets & shps] 


14.0 


Q 1 7Q + 0.021 

U-J- 'U _ 0.017 


1 20 + 

U.iZU „ 0.008 


0.002 


+0.009 
-0.008 


resum 




e+e- [jets & shps] 


22.0 


Q 1 p-1 + 0.015 
U.iOi „ 0.013 


n 1 1 + 0.009 

U.iiS „ 0.008 


0.003 


+0.009 
-0.007 


resum 


m 


e+e- [jets & shps] 


35.0 


U.i^O _ 0.007 


n 1 0Q + 0.008 

U.iZO „ 0.006 


0.002 


+0.008 
-0.005 


resum 


m 


e+e- [(Thad] 


42.4 


0.144 ±0.029 


0.126 ±0.022 


0.022 


0.002 


NNLO 


IHSiiaS] 


e+e- [jets & shps] 


44.0 


Q 1 OQ + 0.011 

U.iOJ _ 0.008 


n 1 90 + 0.008 

U.iZO _ 0.006 


0.003 


+0.007 
-0.005 


resum 




e+e- [jets & shps] 


58.0 


0.132 ±0.008 


0.123 ±0.007 


0.003 


0.007 


resum 


w\ 


pp bbX 


20.0 


n 14c: + 0.018 
U.i^O _ 0.019 


0.113 ±0.011 


+ 0.007 
- 0.006 


+ 0.008 
- 0.009 


NLO 


m 


pp, pp 7X 


24.3 


Q IOC- + 0.012 
U.iOO „ 0.008 


^^^^^^ + 0.008 
U.ilU _ 0.005 


0.004 


+ 0.007 
- 0.003 


NLO 


m 


o-(pP ^ jets) 


40 - 250 




0.118 ±0.012 


+ 0.008 
- 0.010 


+ 0.009 
- 0.008 


NLO 


m 


e+e- r(Z had) 


91.2 


n 1 0.0058 

U.iZZD_ 0.0038 


n 1 00R+ 00058 
U.iZZD_ 0.0038 


±0.0038 


+0.0043 
-0.0005 


NNLO 


m 


e+e- 4-jet rate 


91.2 


0.1176 ±0.0022 


0.1176 ±0.0022 


0.0010 


0.0020 


NLO 


m 


e+e- [jets & shps] 


91.2 


0.121 ±0.006 


0.121 ±0.006 


0.001 


0.006 


resum 


m 


e+e- [jets & shps] 


133 


0.113 ±0.008 


0.120 ±0.007 


0.003 


0.006 


resum 


m 


e+e- [jets & shps] 


161 


0.109 ±0.007 


0.118 ±0.008 


0.005 


0.006 


resum 




e+e- [jets & shps] 


172 


0.104 ±0.007 


0.114 ±0.008 


0.005 


0.006 


resum 




e+e- [jets & shps] 


183 


0.109 ±0.005 


0.121 ±0.006 


0.002 


0.005 


resum 


m 


e+e- [jets & shps] 


189 


0.109 ±0.004 


0.121 ±0.005 


0.001 


0.005 


resum 


m 


e+e- [jets & shps] 


195 


0.109 ±0.005 


0.122 ±0.006 


0.001 


0.006 


resum 


m 


e+e- [jets & shps] 


201 


0.110 ±0.005 


0.124 ±0.006 


0.002 


0.006 


resum 


m 


e+e- [jets & shps] 


206 


0.110 ±0.005 


0.124 ±0.006 


0.001 


0.006 


resum 
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Table 2: Measurements of as(Mzo) included in the process to determine the world average, c.f. ta- 
ble ^ The rightmost two columns give the exclusive mean value of as{Mzo) calculated without that 
particular measurement, and the number of standard deviations between this measurement and the 
respective exclusive mean, treating errors as described in the text. The inclusive average from all listed 
measurements gives as(Mzo) = 0.1189 ± 0.0007. 



Process 


Q [GeV] 




excl. mean as(Mzo) 


std. dev. 


DIS [Bj-SR] 


1.58 


n 1 91 + 0.005 

U.iZi _ 0.009 


0.1189 ±0.0008 


0.3 


r-decays 


1.78 


0.1215 ±0.0012 


0.1176 ±0.0018 


1.8 


DIS [z/; xFs] 


2.8 - 11 


n n q + 0.007 

U.iiJ _ 0.006 


0.1189 ±0.0008 


0.0 


DIS [e/fi; F2] 


2 - 15 


0.1166 ±0.0022 


0.1192 ±0.0008 


1.1 


DIS [e-p — > jets] 


6 - 100 


0.1186 ±0.0051 


0.1190 ±0.0008 


0.1 


T decays 


4.75 


0.118 ±0.006 


0.1190 ±0.0008 


0.2 


QQ states 


7.5 


0.1170 ±0.0012 


0.1200 ±0.0014 


1.6 


e"'"e" [r(Z — > had) 


91.2 


n 1 00fi+ 0-0058 
U.iZZO_ 0.0038 


0.1189 ±0.0008 


0.9 


e"'"e~ 4-jet rate 


91.2 


0.1176 ±0.0022 


0.1191 ±0.0008 


0.6 


e"'"e~ [jets & shps] 


189 


0.121 ±0.005 


0.1188 ±0.0008 


0.4 
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